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Root systems are critical for water and nutrient acquisition by crops. Current methods measuring root 
biomass and length are slow and labour-intensive for studying root responses to environmental stresses in 
the field. Here, we report the development of a method that measures changes in the root DNA 
concentration in soU and detects root responses to drought in controlled environment and field trials. To 
allow comparison of soil DNA concentrations from different wheat genotypes, we also developed a 
procedure for correcting genotypic differences in the copy number of the target DNA sequence. The new 
method eliminates the need for separation of roots from soil and permits large-scale phenotyping of root 
responses to drought or other environmental and disease stresses in the field. 

Drought (water deficit) represents the most severe constraint to agricultural production, and accounts for 
approximately 70% potential yield loss worldwide''^. Root systems are critical for uptake of water and 
nutrients by plants. Effective deployment of genetic variation in root response to drought could have 
significant implications for world food securityl 

Studies on growth of root systems in the field have been restricted by lack of high throughput methods'*. 
Current methods for high throughput phenotyping are based on measurements of seedlings grown in artificial 
media such as agar"""'. However, root growth in artificial media can differ from that of plants grown in hetero- 
geneous soil environments'*''' Characterising root growth in the field is challenging. Most studies measure 
changes in root biomass and length, which requires roots to be extracted from soil, separated from the dead roots 
and roots of other species, and then quantified using image analysis, and/or dried to measure root dry weighf*. 
Other methods such as excavation and minirhizotron have also been used" '^. 

Cereals such as wheat have fibrous root systems, and the fine lateral roots (approximately 0.07-0.15 mm in 
diameter) account for most of the total root length*''. These fine roots fall below the current detection limit for X- 
ray computed tomography imaging*'*. They are usually fragile and translucent, and difficult to recover for accurate 
measurements*^ *^, which limits their application to phenotype root growth for field-grown plants'*. 

A decade ago specific plant DNA sequences were amplified using polymerase chain reaction (PCR)*" for 
qualitative identification of mixed tree roots isolated from soils*'. Recently, quantitative real-time PCR (qPCR) 
has been used to measure plant DNA concentrations from isolated roots of wild grasses of single or mixed species 
for estimation of relative abundance of roots of pot-grown plants™"^^. qPCR has also been used to measure plant 
DNA concentrations directly from soil for estimation of root dry weight of different forage grasses and root 
distribution in soil profiles for field-grown plants^^"^'*. 

Plant DNA in soil samples is largely derived from live roots. The remnant DNA in dead roots is expected to 
contribute little to the total DNA concentration as over 90% of plant DNA in soils degrades within 7 to 10 days 
after plant death^^'^**. Plant DNA concentrations in soils can provide a measure of root response to treatments such 
as addition of lime to acid soU^^ and complete defoliation or a herbicide treatment^^, but also of root distribution in 
soil profiles^*'^'*. 

The measurements of plant DNA concentration in soils with qPCR are through amplification of a target DNA 
sequence. The internal transcribed spacers {ITSl and ITS2) of the nuclear ribosomal DNA (rDNA) unit are widely 
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used for phylogenetic studies and DNA barcoding as they show high 
levels of interspecific divergence^^'^''. ITS2 sequences have been used 
as a target sequence for qPCR measurements of plant DNA concen- 
tration of different forage grass species in soiP^'^* because of its high 
specificity at the species level and high copy numbers^' allowing 
detection of low DNA concentrations^^. There have been no reports 
so far of genotypic differences in root response to environmental 
stresses and root distribution within crop species, using plant DNA 
concentrations. In this report, we describe the development of a 
method for quantitative measurement of wheat DNA concentrations 
directly in soils. Since there have been reports of interspecies variation 
in ITS copy numbers^^ and intraspecies variation in ITS sequences^*"'^", 
we devised a procedure for correcting genotype differences in the 
copy number of the target DNA sequences. We demonstrate that 
the DNA-based method provides a new tool for high-throughput 
measurements of root responses to drought in field trials. 

Results 

Genotypic variation in copy numbers of ITS2 sequences. Genoty- 
pic differences in copy numbers of wheat target sequence ITS2 
(TalTSl) among 20 modern Australian spring wheat genotypes 
were determined using qPCR of TalTSl and TaPHTl;6A, a single- 
copy gene, (Supplementary Fig. 1). Ct values of TaITS2 and TaPH- 
T1;6A from root DNA samples were used to calculate TaITS2 copy 
number per genome. The Ct values of TaPHTl;6A were similar 
among 20 genotypes (27.68 ± 0.06), whereas the Ct values of 
TaITS2 among 20 genotypes with sequence-specific primers, 
ranged from 13.82 ± 0.20 to 16.16 ± 0.18. The genotypic diffe- 
rences in the Ct value of TaITS2 were highly significant (P < 
0.001). Calculated copy numbers of TaITS2 varied from approxi- 
mately 2,900 copies per genome in the variety Seri to 15,000 copies 
per genome in Young (Fig. 1). The average copy number of ra7TS2 in 
the 20 genotypes was approximately 8800 per genome, comparable 
to 8200 rDNA repeat units per genome estimated for the two major 
rDNA loci (Norl and Nor2) in wheat by DNA hybridization^'. The 
results indicate that there is genotypic variation in TaITS2 copy 
numbers among the 20 wheat genotypes. 



As other plant species closely related to wheat {Triticum aestivum) 
could be present in the field, such as T. urartu, T. monococcum, Aegi- 
lops speltoides, A. tauschii and Secale cereal, amplification refractory 
mutation system (ARMS)"*" primers were designed for TaITS2 with 
deliberate 3' mis-matches to eliminate DNA amplification from 
these related plant species. The ARMS primers were tested for effi- 
cacy in assessment of genotypic variation in TaITS2. The mismatches 
in the ARMS primers increased Ct values of TaITS2 for all 20 wheat 
genotypes, but there was a highly significant correlation (r^ = 
0.986) in copy numbers of TaITS2 per genome between ARMS 
and sequence-specific primers (Supplementary Fig. 2). This result 
further confirms the genotypic variation in TaITS2 copy numbers 
among the 20 wheat genotypes, and indicates that the ARMS primers 
are suitable for amplification of TaITS2. The ARMS primers were 
therefore used in DNA quantification for all soil samples. 

For genotype comparison, it was necessary to correct the differ- 
ences in copy numbers of TaITS2 among 20 genotypes used in the 
experiments. The copy number of TaITS2 in Krichauff was selected 
as the reference genotype because it was close to the average copy 
number of TaITS2 per genome (Supplementary Table 1). A scaling 
factor was calculated for each genotype to correct the difference in 
the copy number of TaITS2 (Supplementary Table 1). 

Genotypic variation in root responses to drought in controlled 
and field environments. Root responses to drought were initially 
studied in pot-grown plants in a controlled environment. As root 
length density (RLD) is often used to measure root growth in pot and 
field-grown plants^', both RLD and root DNA density (RDD) were 
examined in two wheat cultivars. Kukri and Gladius for root respon- 
se to stable and moderate drought ( — 0.4 MPa soil water potential). 
The drought treatment started after seedling establishment and 
ended at flowering. Genotypic differences in shoot dry weight at 
flowering were significant in the well-watered plants between 
Kukri and Gladius, but not in the drought-treated plants (Fig. 2a). 
The drought treatment reduced shoot growth by 35% in Gladius and 
53% for Kukri relative to the well watered controls (Fig. 2a). Diffe- 
rences in RLD were similar to shoot dry weight for the genotypes and 
water treatments (Fig. 2b). Kukri had significantly higher RLD than 
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Figure 1 | Genetic variation in copy numbers of TalTSZ per genome among 20 wheat genotypes. Total root DNA was isolated from primary roots of one 
month-old plants grown in a potting mix, and used to determined copy numbers of TalTSl per genome using quantitative real-time PGR. Ct values of 
TaITS2 and TaPHTl;6A using sequence-specific primers and probes were determined for all samples. The copy number of TaITS2 per genome was 
calculated by 2*"^' (ACt = TaPHTl;6A Ct - TalTSl Ct). Predicted means of four biological replicates are presented. 
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Figure 2 | Effect of stable and moderate drought on shoot dry weight, root length density and root DNA density on wheat plants grown in pots under a 
controlled environment. Effect of stable and moderate drought on shoot dry weight (a), root length density (b) and root DNA density (c) of two wheat 
cultivars, Kukri and Gladius. Soil moisture tension in the stable and moderate treatment was maintained at — 0.4 MPa for four weeks after sowing until 
flowering. Well watered (WW) plants were the control. Plants were grown in a controlled climate environment. Means of five replicates are represented. 
There are significant differences in shoot dry weight and RLD for interactions of genotype X drought {P < 0.05). The differences in RDD for the 
genotypes and drought treatments were significant (P < 0.006), but the interaction of genotype X drought was not significant (P = 0.56). Vertical lines 
indicate the least significant difference at P = 0.05. 



Gladius in the well-watered plants, but no difference in RLD was 
observed between the two genotypes exposed to the drought treat- 
ment (Fig. 2b). The reduction in RLD by the drought treatment was 
46% in Gladius and 60% in Kukri (Fig. 2b). In contrast, Gladius had 
significantly higher overall RDD than Kukri regardless of the water 
treatments (Fig. 2c). The drought treatment significantly reduced 
RDD by approximately 30% in both genotypes relative to that of 
the well watered treatment (Fig. 2c). These results indicate that 
RDD detects root responses to the drought as RLD, and reveals the 
difference between RDD and RLD. 

Next, root responses to a Mediterranean-type drought scenario""^ 
was measured using RDD. Plants of the two wheat cultivars. Kukri 
and Gladius were grown in field soil under a rain-out shelter, and 
exposed to two cycles of drought during flowering (Fig. 3a). The 
cyclic drought reduced soil moisture potential to —0.5 MPa 
(Fig. 3a), leading to a significant reduction of RDD in both cultivars 
(at least 64% relatively compared to well watered plants) in the top 
10-cm soil layer (Fig. 3b). There was no effect of the cyclic drought on 
RDD of either cultivar below 10-cm depth (Fig. 3b). This indicates 
that RDD in the top 10-cm soil is particularly sensitive to the cyclic 
drought. RDD of both genotypes in the well watered and cyclic 
drought plots was significantly lower in 20-30 cm layer than in the 
top soil layers (Fig. 3b). Furthermore, RDD of Gladius across the 
three depths and the drought treatments was significantly higher 
than that of Kukri, revealing the genotypic variation in response to 
the cyclic drought between the two genotypes. 

The vertical distribution of roots from 0 to 40 cm depth was also 
measured using RDD at flowering for two wheat cultivars, Krichauff 
and Berkut during a severe drought in 2007. RDD of the two cultivars 
was the greatest in the 0-10 cm layer (Fig. 3c), similar to the well- 
watered plants (Fig. 3b). There was a sharp decline in RDD of both 
cultivars below 10-cm depth, and little change from 10 to 40 cm 
depths (Fig. 3c). However, Berkut had a significantly higher RDD 
than Krichauff across all four depth measurements (Fig. 3c). These 
results reveal that there is the genotypic difference in vertical distri- 
bution of live roots under a severe drought condition. 

High-throughput screen for genetic variation in root response to 
drought. To test whether the DNA-based method is suitable for 
studying genetic variation in root response to drought, a large trial 
consisting of twenty modern wheat genotypes with diverse genetic 
background was conducted at five locations in South Australia in 
2008. The growing- season rainfall in 2008 was well below the yearly 
average and each location varied in accumulated rainfall of the 



growing season (Supplementary Table 2). On average, there was an 
approximately two-fold variation of RDD in the top 10-cm soil at 
flowering across the 20 genotypes for the five locations (Fig. 4a). 
Genetic variance at each location was highly significant, and the 
heritability at each site was greater than 0.5 (Supplementary Table 
3), indicating a strong genetic component underlying RDD. The 
genetic correlations between sites were also highly significant 
(Supplementary Table 3). 

Genotypic variation in RDD was also examined for root vertical 
distribution at flowering at one site, Karoonda. Seven of the 20 gen- 
otypes were assessed to 45-cm depth (Fig. 4b). No significant differ- 
ences in RDD across three depths were observed for Gladius and 
Berkut, whereas significant differences were seen between the 0- 
15 cm layer and those deeper in the soil profile for the other five 
cultivars (Fig. 4b). These results indicate that different wheat geno- 
types have distinct root responses to severe drought in terms of RDD 
distribution in the soil profile. 

Discussion 

Current methods based on root biomass and length can only provide 
low-throughput measurements for root growth in field-grown crops, 
which restricts studies on genetic diversity of root systems in the 
field**. Quantitative PGR allows measurements of plant DNA con- 
centration in soil and can detect short-term and long-term changes 
of root systems in responses to stresses such as drought (Fig. 2c, 
3b, 3c). Overall, these data indicate that plant DNA concentration 
in soil can be used to assess root responses to stresses, which may not 
be obtainable from root biomass or length. Plant DNA concentra- 
tions in soils can also be used to detect wheat root distribution in soil 
profiles (Fig. 3c, 4c). The vertical distribution of RDD are generally 
analogous to that measured with root biomass or length in field- 
grown wheat plants'"" ''*. In addition, the plant DNA concentration 
in soil can be used for estimation of root biomass and length with 
careful calibration^'"^'*. However, the proportion of live root cells in 
soil changes with different field sites, soil depths, and ages of plants^'. 
The calibration is needed for each situation as it is unlikely that there 
exists a close correlation between plant DNA concentration and root 
biomass or length in different soils. 

Two types of target sequences (inter simple sequence repeat and 
ITS) have been used to quantify plant DNA concentration by 
qPCR^""^'. ITS2 sequences were used in this study as they are highly 
specific at the species level and are in high copy numbers^^, which 
allows detection of low DNA concentrations in soiP^. However, 
potential differences in the copy number of the ITS2 sequences 
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Figure 3 | Effects of cyclic drought and continuous drought on root DNA 
density (RDD) in soil profile, (a) Soil moisture tension was measured at 
15 cm and 35 cm depth. Two cycles of the drought treatment (Drought) 
were applied to the plants grown in a rain-out shelter at 55 and 65 days after 
seed sowing, respectively. Well watered (WW) plants were the control, 
(b) Effects of the cyclic drought on RDD of two wheat genotypes, Kukri 
and Gladius in 0-30 cm depth. Four soil cores (2.5 cm in diameter) were 
collected at 70 days after seed sowing from between rows for each plot. 
Means and standard errors of four replicates are presented for RDD. There 
are highly significant differences in RDD for cultivars, cyclic drought and 
depth, and for interactions of cyclic drought X depth X genotypes (P < 
0.003). (c) RDD of two wheat genotypes, Krichauff and Berkut in 0-40 cm 
depth at flowering. Five soil cores (2.5 cm in diameter) were collected from 
between rows of plots at Roseworthy, South Australia in 2007. Means and 
standard deviations of three replicates are presented for RDD. There are 
highly significant differences in RDD for cultivars and depth (P < 0.001). 
Ct values of TaITS2 were determined for all soil samples using quantitative 
real-time PGR with ARMS primers and the TaITS2 probe. RDD presented 
was standardised with the scaling factor for each genotype. 



between genotypes need to be quantified and corrected for genotype 
comparison. Such correction allows comparison in plant DNA con- 
centration in soU between wheat genotypes (Fig. 2c, 3b, 3c, 4a, 4b). 
There was a strong genetic correlation in RDD across multiple field 
sites (Supplementary Table 3) indicating that RDD is a heritable trait. 
In addition, the results from the field trials demonstrate that the 
DNA-based method can provide high-throughput assessment of 
root responses to drought in the field. These results show that 
RDD could assist in the characterisation of genetic diversity in root 
system architecture and root responses to water and nutrient supply 
in field-grown crops. 

The association of RDD with drought tolerance in yield has 
been observed for some of the genotypes (Fig. 4a). The southern 
Australian or Mediterranean-type environment is characterised by 
late season drought and low soil moisture and poor quality. In this 
environment deep rooting offers few advantages since there is little 
useful moisture at depth. In much of southern Australia, ground- 
water is saline or there are toxic levels of boron at depth. Genotypes 
that perform well in these environments'^, Excalibur, RAC875, Yitpi 
and Wyalkatchem had higher RDD in the top 10 cm soil, while less 
well adapted genotypes, such as Correll, Drysdale, Kukri and Young 
had lower RDD. However, two drought tolerance cultivars, Gladius 
and Mace show relatively low RDD (Fig. 4a). It may be premature to 
draw definitive correlation between RDD and yield. Yield is deter- 
mined by a complex interaction of a wide range of environmental 
parameters over the full developmental phases of plant growth. 
Consequently, further analyses are required to separate effect of root 
traits in terms of RDD on yield from other traits. As more data are 
accumulated on changes in RDD across different germplasm pools 
and different environments, the significance of RDD relative to other 
measures of stress tolerance should become clearer. In addition, 
future studies should compare RDD with results from X-ray com- 
puted tomography" and other in situ techniques for live roots, such 
as a green fluorescent protein reporter'". 

In summary, the measurement of wheat root DNA in soil by qPCR 
can detect root responses to drought with high throughput. It is 
expected that this method will provide a new tool to complement 
other techniques to characterise root responses to abiotic and biotic 
stresses in field conditions. The application of this new technique will 
promote exploration of genetic diversity in root systems for improve- 
ment of water and nutrient acquisition''''"'. This new technique 
should be readily applicable to crops other than wheat and diverse 
field environments and natural ecosystems. 

Methods 

Determination of genetic variation in TaITS2 copy numbers among 20 wheat 
genotypes. Twenty wheat {Triticum aestivum L.) genotypes were grown in pots 
(15 cm in diameter X 15 cm in height) which contained 2 kg of University of 
California potting mix with complete plant nutrients. One pre-germinated seed per 
pot was grown in a growth room with 15/10"'C (day/night) temperature and 12 hr 
day/night cycle for one month. Light density was 300 [rmol m"** s"' at the plant level. 
At harvest, pots were soaked briefly in water, and the entire root system was carefully 
washed free of potting mix. Excess water on roots was blotted on laboratory tissues. 
Primary roots with laterals were collected, immediately frozen in liquid nitrogen, and 
freeze-dried. 

Total DNA from freeze-dried roots was isolated as described for soil DNA 
extraction (see below). DNA concentrations of samples were determined using 
QuantiFluor™ dsDNA System (Promega, Madison, WI, USA), and diluted to 
2.5 ng fiL"\ Quantitative real-time PGR using TaqMan® Minor Groove Binder 
(MGB) probes (Applied Biosystems, Foster Gity, GA, USA) was performed on ABI 
PRISM® 7900HT as described'*'^. For TaITS2, the Gt value was determined using the 
TaqMan probe (6FAM 5'-CGGGATGTGGTGGGT-3') with two different primer 
pairs. TaITS2 sequence- specific primer pair consisted of forward primer (5'-CGAA- 
AAGAGGCTGCGAACGA-3') and reverse primer (5'-TGGGTCGAGCGTGC- 
CTTGA-3'). The AI^S primer pair for TalTS2 consisted of forward primer (5'- 
CGAAAAGAGGGTGGGAAACA-3') and reverse primer (5'-TGGGTGGACGGT- 
CGCTAGA-3' ). The ARMS primers were designed with 3' mis-matches for increased 
specificity to discriminate T. aestivum closely- related species. For single copy gene, 
TaPHTl;6A located on chromosome SAL (Supplementary Fig. 1), Gt values were 
determined using forward primer (5'-CGGACGTGGTCGGAGAGT-3'), and reverse 
primer (5'-AAAGCATGAAAAGGATCGATAGAC-3') with the TaqMan® MGB 
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Figure 4 | Genetic variation of root DNA density (RDD) in field-grown wheat genotypes, (a) Genetic variation of RDD in top 10-cm soil among 20 
wheat genotypes at five field sites. Ten representative soil cores (one cm in diameter X 10 cm in depth) from between rows were collected at flowering for 
each plot. Predicted means of four replicates are presented for RDD. (b) Genotypic variation of RRD in soil profile among seven wheat genotypes. 
Five representative soil cores (2.5 cm in diameter and 45 cm in depth) were collected at flowering from between rows for each plot at one field site, 
Karoonda. The 45-cm soil cores were separated into three 15-cm sections (0-15 cm, 15-30 cm and 30-45 cm), and the corresponding sections of five soil 
cores were combined for root DNA analysis. Natural logarithm transformed data were used for statistical analysis and means (n = 3) are presented for 
RDD. There are significant differences in RDD for genotype and depth (P < 0.001 ), and for interactions of genotype X depth (P < 0.034). Vertical lines 
indicate the least significant difference (LSD0.05) for interactions of genotype X depth. Ct values of TaITS2 were determined for all soil samples using 
quantitative real-time PGR with ARMS primers and the TaITS2 probe. RDD was standardised with the scaling factor for each genotype. 



probe (6FAM 5'-CCAACACCTGCATAGTA-3' on the reverse strand). For each 
wheat genotype, the copy number of TaITS2 per genome was calculated by 2^*^' (ACt 
- TaPHTl;6A Ct - TaITS2 Ct). 

Measurements of root DNA density in soil. Soil cores were dried overnight at 45 ''C 
in a dehydration oven, and DNA was extracted by the commercial service. Root 
Disease Testing Service (SARDI, Adelaide, SA, Australia). This DNA extraction 
method shows consistency with two other methods available for root DNA isolation 
and measurements^^. Quantitative real-time PCR was used to determine Ct values of 
soil samples using the ARMS primer pair of TaITS2 with the TaqMan probe of 
TaITS2 as described above^^. The Ct values could not be directly compared in 
different batches of samples. Therefore, the Ct values were converted to pg root DNA 
g~* soil using standard calibration curves obtained from 10-fold serial dilution of 
Krichauff root DNA, ranging from 200 to 0.002 pg jil"' for TaITS2 and from 20,000, 
to 2.0 pg for TaPHTl;6A measurements. The root DNA of Krichauff was isolated 
from 7 day-old, sterile roots. 

Plant accelerator experiment. Two genotypes (Gladius and Kukri) were used in the 
experiment. Plants were grown in white plastic pots (14 cm in diameter X 19 cm in 
height) containing 2.5 kg of soil mixture (equal volume of clay-loam soil and coco- 
peat mix) with complete fertilizers. Three seeds per pot were sown on 2P' February, 
2011, and thinned to one plant per pot two weeks after sowing. The plants were 



treated with two watering regimes from 28 days after sowing. In the well-watered 
regime, soU water content was maintained at 20% by watering every second day; 
equivalent to a soil water potential of approximately — 0.05 MPa. Soil water potential 
was estimated by weighing pots every two days. For the moderate drought regime, 
watering was stopped at 28 days after sowing. When soil water potential reached to 
— 0.4 MPa, the pots were rewatered. This mean level of soU water potential was 
maintained untU the end of the experiment by watering every second day. 

Plants were grown in the glasshouse (phenotyping platform of "The Plant 
Accelerator", Adelaide, Australia). Air temperature and relative humidity were 18^ 
day/24 "C night, and 50-90% humidity. Plants were harvested at 52 days after sowing, 
24 days after being exposed to the two water regimes. Shoot dry weight was recorded 
after 3 days at 85"C. All soU in the pot was carefully removed and split into four 
quarters, two for analysis of RDD and the other two for root length density. 

Root DNA density in a rain-out shelter experiment. Plants of two wheat cultivars. 
Kukri and Gladius were grown in small experimental plots (20 X 60 cm, 8 plants per 
row, two rows per cultivar) in a rain-out shelter at the Waite Campus, Adelaide, 
Australia. The plots were fertilized twice with 8 g Aquasol m~^ (Hortico, Clayton, 
VIC, Australia) at 37 and 50 days after sowing. Seeds were sown on 3rd September, 
2009. Two cycles of drought treatments were applied to plants at 55 and 65 days after 
sowing. Two moisture sensors were embedded at 15-cm and 35-cm depth, 
respectively to monitor soU moisture potential. Four soU cores (30 cm long X 2.5 cm 
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wide) were taken at 70 days after sowing from between rows in each plot for each 
wheat cuhivar. The 30-cm soil cores were divided into three 10-cm sections, and the 
corresponding sections of four soil cores from each plot were combined for analysis of 
root DNA density. 

Root DNA density in field experiments. Two field experiments were conducted. The 
first was at Roseworthy {South Australia) in 2007. The distribution of plant DNA in 
the soil profiles was examined for two wheat mapping parents, Berkut and Krichauff. 
Five cores (40 cm long X 2.5 cm wide) were taken from between rows of three plots 
per cultivar, and then divided into 10-cm sections for analysis of root DNA density. 
The second multiple genotype experiment was conducted at five sites in South 
Australia in 2008. The growing- seasonal rainfall in 2008 was the lowest in decades. 
Rainfall data were obtained from weather stations closest to the trial sites 
(Supplementary Table 2). Twenty-one, hexaploid semi-dwarf wheat genotypes from 
diverse genetic background and widely used in breeding programs around Australia 
were included. The trial design was randomized split block design (12 ranges X 7 
rows), and four replicates. Plots were 3.20 X 1.25 m, row spacing 0.25 m, and 
managed as per district practice. Seeding rate was 200 seeds per square metre. Soil 
cores (10 cm long X one cm wide) were collected from the top soil at flowering. 
Twenty cores were taken from between rows, and combined to make a composite 
sample per plot for plant DNA analysis. As mixed seeds were found for one of 21 
genotypes, RDD data of 20 genotypes were presented. 

For the genotypic variation in the distribution of roots in the soil profile, soil cores 
were taken at Karoonda site from between rows of seven wheat genotypes with three 
replicates. Five soil cores (45 cm long X 2.5 cm wide) per plot were collected at 
flowering four days after the top soil sampling. The 45-cm cores were separated into 
three 15-cm sections, and the corresponding sections of five cores were combined to 
make a composite sample for plant DNA analysis. 

Root length measurements. Soil cores were soaked in water. Roots were carefully 
washed free of soil particles, and plant debris was removed from roots. The roots were 
spread out in root position tray (Reagent Instruments Inc., Canada) and scanned at 
600 dpi. Scanning parameters were set as those described^^. Length of roots was 
determined using the WinRHIZOPro Image Analysis software (Reagent Instruments 
Inc., Canada). 

Statistical analyses. Experiments were analysed using the R (R Development Core 
Team, 2011, URL http://www.R-project.org/) package ASReml-R^^. In the statistical 
analyses of TaITS2 copy number per genome, shoot dry weight, RLD and RDD from 
the pot experiment in controlled environment, RDD from the rain-out shelter 
experiment, the first field experiment, and deep soil cores in the second field 
experiment, genotypes were fitted as fixed effects, along with the other treatments of 
interest, with an analysis of variance used to test the significance of effects. For RDD of 
the 0-10 cm soil cores from the second field experiment, a linear mixed model was 
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